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It  is  shown  that  cubic  spline  interpolation  with  the  not-a-knot  side 

y 

condition  converges  to  any  C  -interpoland  without  any  mesh-ratio  restriction 

/ 

as  the  mesh  size  goes  to  zero.  \ 
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SIGNIFICANCE  AND  EXPLANATION 


The  cubic  spline  interpolant  to  a  given  function  agrees  with  that 
function  at  the  knots  t^  <  ...  <  tR  of  the  splines  its  complete 
specification  requires  two  additional  conditions.  Originally,  the  slope  at  the 
end  points  t1  and  tR  was  used,  giving  the  cas^lete  spline  interpolant. 

But,  in  the  absence  of  such  information,  one  has  to  try  something  else.  One 
proposal  is  the  not-a-knot  condition  [B3;Ch.  IV]  in  which  the  third  derivative 
of  Pg  is  made  continuous  at  t2  and  at  tn_1  .  In  effect,  these  next-to- 
boundary  knots  are  not  knots.  But,  whereas  the  convergence  of  complete  spline 
interpolation  has  been  established  in  various  ways,  there  doesn't  seem  to  be  a 
proof  of  the  convergence  of  not-a-knot  cubic  spline  interpolation  in  the 
literature. 

The  present  report  remedies  this  situation. 
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CONVERGENCE  OF  CUBIC  SPLINE  INTERPOLATION 


WITH  THE  NOT-A-KNOT  CONDITION 
Carl  da  Boor 

Tha  cubic  aplina  intarpolant  Pg  to  a  given  function  g  agrees  with  g  at  the 
knots  tj  <  ...  <  tn  of  tha  spline;  Its  ccoplete  specification  requires  two  additional 
conditions.  Originally,  the  slope  at  tha  end  points  t^  and  tn  was  used,  giving  the 
collate  aplina  intarpolant.  But,  in  the  absence  of  such  information,  one  has  to  try 
soaiathing  else.  One  proposal  ia  the  not-s-kaot  condition  [B3;Ch.  IV]  in  which  the  third 
derivative  of  Pg  is  made  continuous  at  t2  and  at  tn_j  •  In  effect,  these  next-to- 
boundary  knots  are  not  knots.  But,  whereas  the  convergence  of  complete  spline 
interpolation  has  been  established  in  various  ways,  there  doesn't  seem  to  be  a  proof  of 
the  convergence  of  not-a-knot  cubic  spline  interpolation  in  the  literature  (except  for 
[SV]  in  which  it  appears  as  a  special  case  but  is  only  treated  for  a  uniform  mesh). 

The  argument  here  follows  the  standard  path  (see  [B,],  [B21):  One  considers  the 
derived  projector  P"  obtained  from  P  by  the  rule 

(1)  P"D2g  D2(Pg) 

and  shows  it  to  be  bounded  in  the  uniform  norm  independently  of  the  meBh,  i.e.,  one  shows 
that 

(2)  N  supt  «P"I  <  - 
with 

IP-1  supf  lP"fl/lf« 

and 

If!  i“  I f  *  —  sup  |  f ( x ) |  . 

t  <x<t 
1  n 

Since  P"  is  a  linear  projector  with  range 


i.e,  with  range  the  span  of  the  linear  B-splines  on  the  knot  sequence 

3  <  si  >°+2  :»  . tn-3'tn-2'tn'tn>  *  henc«  m  n“2  ' 

this  implies  that 

lD2(g  -  Pg)  I  =*  lD2g  -  P"(D2g)l  <  I1-P"l  dist(D2g,S)  <  (1+M)  a  (|s|) 

D  g 

with  the  modulus  of  continuity  of  f  (on  [tf,tn)>  and 

| s  |  s«  rnaxj^  As^  . 

Consequently, 

lg  -  Pgl  <  t(1+M)|t|2/8}(|s|2/8)lD4gl  , 

using  the  fact  that  Pg  agrees  with  g  at  t  and  using  the  standard  estimate  for 

2  2  2 

dist(D  g,S)  in  case  D  (D  g)  is  bounded. 

This  leaves  the  hard  part  of  the  argument,  namely  a  proof  of  (2),  but  here,  too,  the 
approach  is  standard,  as  follows. 

2  2 

Since  Pg  agrees  with  g  at  t  ,  P"D  g  agrees  with  D  g  at 

Xi  !=  ^  Mi*  '  1=1 . m  ' 

with  :=  j  t  the  linear  B-splines  for  the  knot  sequence  t  (and  normalized  to 

have  j  «  1  ).  This  follows  from  the  fact 

2!  [ti,ti+,,ti+2]g  *  /tMx)D2g(x)dx  . 

Hence,  with  :=  Ni  2  s  the  linear  B-splines  for  the  knot  sequence  s  (and  normalized 

to  sum  to  1  ),  and  f  :«  D2g  ,  we  obtain  P"f  in  the  form 

P"f  =:  Z  at  , 

where  a  is  the  unique  solution  of  the  linear  system 

(3)  As  =■  ( *>if ) 

with 

(4)  A  :=  t  )  . 

Consequently, 

«P"f«  =  '»•„  <  IA_1I  jnaxi|/Mlf  |  <  «A-1l<.lfl<.  . 

This  shows  that 

(5)  <P"«  -  'A_1»„ 

and  so  reduces  the  proof  of  (2)  to  finding  a  uniform  lower  bound  for  A  . 
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Since  A  is  totally  positive,  a  lower  bound  for  A  is  obtained  fron  any  vector  x 
for  which 

d  min j  (-jllAx)^  >  0  . 

For,  the  total  positivity  of  A  implies  (by  Cramer's  rule)  that  its  inverse  is 
checkerboard,  i.e., 

(-)i+iA_1(i,j)  >  0  . 

Therefore, 

|x(i)|  -  | £A“1(i, j) <Xx>  ^ |  -  I  |A"1(i,j)|(-)j(Ax)j  >  l|A-1(i,j)|  d 

hence 

•  xl,.  ■  maxi|xi|  *  lA-1!^ 

or 

(6)  'A*’11..  <  • 

The  matrix  3A  is  tridiagonal,  with  typical  row 

Asi/<8i+2"“i>  ,  2  ,  ^s^+ ^  '  (“3, •* » ,#“2  , 

and  only  the  first  two  and  the  last  two  rows  deviate  from  this.  Here  are  the  first  two 
rows: 


1  +  (At^/ls.j-S,) 


2  -  (AtjJ/tS^S,) 


(At,) 


2<at„r  +  34t,dt, 


(s. 


^vv  ■  lt;vi  +  -Si  ' 


The  last  two  rows  can  be  obtained  from  this  by  the  appropriate  change  in  variables. 

The  typical  row  of  3A  suggests  the  vector  x  : “  (“1,1, -1,1,  ...  )  as  suitable  since 
this  gives 

b  :■  3 Ax  *  ( b^ , b2 , ”1 , 1 ,  —  1 , 1 , • • • ,b^_ ^ , b^)  , 

i.e.,  an  appropriate  result  except,  perhaps,  for  the  first  two  and  last  two  components. 

One  can  check  that,  as  t2  varies  in  the  interval  [s^Sj],  the  first  component  of  b 
becomes  positive  and  this  destroys  the  desired  alternation.  This  can  be  helped,  though,  by 


choosing  x.  and  x,  in  dependence  on  t2 


Specifically,  X  choose  x  in  the  form 


■o  that 


b  ■  ...  )  i 

with  tha  analogous  happaning  at  tha  othar  and.  This  requires  that  m  >  5  ,  to  avoid 
destructive  interference  of  the  machinations  on  one  end  with  those  at  the  other.  With  this 
assumption,  though,  it  ie  sufficient  to  run  through  the  details  for  the  left  end  only. 

This  means  that  I  solve  the  linear  system 


A(  1, 1  )x j  ♦  A(1,2)k2  -  -1 

A(2,1)x1  +  A(2,2)X2  -  1  -  A(2,3)*(-1) 


Since  the  system  is  invariant  under  a  linear  change  in  the  independent  variable,  I 
normalize  the  situation  to  one  where 


s,  -  0,  a3  -  1 

and  define 


k  i«  *t2  «  h  s“  • 

This  means  that  I  am  interested  in  the  range  k  f  (0,1),  h  >  0  .  I*  these  terms,  the 
system  reads 


1+k  2-k 

k2  3k-k2+2h 


<-1,k  ♦  2h )  . 


Its  solution  is 

x  »  (-5k  +  2k2  -  6h  ♦  2hk,  k  +  2k2  +  2h  ♦  2hk)/D  , 


with 


D  (3k  +  2h  ♦  2hk>  *  0  . 

Now  observe  that,  for  (k,h)  e  (0,1]*^  , 

x^_  (-3)  .  (4k  +  2k2  +  8hk)/D  >  0  with  equality  only  if  k  ■  0  , 
-1  -  x1  »  2( 1-k) (k  ♦  2h)/D  *  0  with  equality  only  if  k  »  1  . 

Finally, 

x,  -  1  -  2k(1-k)/D  <  1  ,  with  equality  only  if  k  »  0  or  1  . 


I  conclude  that  x,  e  (-3,-1)  and  x2  e  [2/3,1]  ,  hence 

lxl„  <  3  . 

Also,  bj  <  -1  ,  since  x2  *  1  ,  while  -  (-)*  otherwise,  hence 

3d  isin^  «-)i3bi  >  1  . 

This  shows  that 

(7)  'A-1*.  <  9 

independently  of  t  (as  long  as  n  >  6  ). 

Remarks,  (i)  As  t2  — >  tj  and  tn_^  — >  tn  ,  not-a-knot  cubic  spline  interpolation 
reduces  to  complete  spline  Interpolation,  hence  P"  becomes  least-squares  approximation 
from  $2  t  .  The  above  argument  then  reduces  to  the  standard  one.  In  particular,  the 
computed  x  is  just  the  alternating  vector  (-1,1,-1,1,...  )  and  the  resulting  bound  is 
the  customary  onei  li  «  3  . 

(ii)  It  may  be  possible  to  carry  out  the  argument  by  perturbation,  starting  off  with 
the  known  stability  of  P“  for  complete  cubic  spline  interpolation  and  showing  that  the 
change  in  the  aide  conditions  to  the  not-a-knot  conditions  is  gentle  enough  (at  least  for 
large  n  )  to  change  <P" •  by  a  bounded  amount. 

< lit )  When  I  wrote  this  note  in  August  of  1984,  I  wrote  that  "It  should  be  possible 
to  show  the  boundedness,  independently  of  t  ,  of  the  derived  linear  projector  P'  ,  given 
by  the  rule 

P'(Dg)  D(Pg)  . 

But  the  argument  via  total  positivity  would  be  a  bit  messier."  In  the  meantime,  R.  K. 

Beatson,  who  had  earlier  brought  my  attention  to  the  fact  that  there  did  not  seem  to  be  a 

convergence  proof  for  spline  Interpolation  with  the  not-a-knot  condition,  has 

independently  studied  this  question  and  shown  that,  strictly  speaking,  P'  is  not  boundable 

independent  of  the  mesh.  The  difficulty  lies  in  the  first  and  last  data  interval.  If  one 

measures  the  size  of  P'f  by  the  number 

IP'fl*  sup  |P’f(t)|  , 

2  n-1 
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